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Abstract 

Identified particle elliptic flow results are presented for the Au+Au reaction at ^snn = 200 GeV as a 
function of transverse momentum and pseudorapidity. Data at pseudorapidities jj ~ 0, 1, and 3.4 were 
obtained using the two BRAHMS spectrometers. Differential V2{rj,pt) values for a given particle type are 
found to be essentially constant over the covered pseudorapidity range, in contrast to the integral V2 values 
which have previously been observed to decrease at forward rapidities. A softening of the particle spectra 
at forward angles is found to account for at least part of the integral V2 falloff. The data are found to be 
consistent with existing constituent quark scaling systematics. 
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1 Introduction 

The azimuthal anisotropy of particle production 
in relativistic heavy-ion reactions is believed to re- 
flect pressure gradients of the thermalized medium 
formed early in the collisions. In a Fourier expansion 
of the azimuthal behavior, the 2 nd -harmonic term, 
with coefficient i> 2 , defines the "elliptic flow" of the re- 
action. Measurements near mid-rapidity at RHIC en- 
ergies have shown that elliptic flow for particles below 

2 GeV is close to that expected from hydrodynamic 
models that assume the formation of a strongly inter- 
acting quark-gluon plasma^ 1 '— . However, integral 
i! 2 values, which are based solely on the asymmetry 
of the overall charged-particle production, are found 
to fall off going away from mid-rapidity, dropping by 
about 35% by pseudorapidity i] = 3^ . 

BRAHMS has now measured i> 2 and the associ- 
ated particle spectra as a function of transverse mo- 
mentum for pions, kaons, and protons at angles (pseu- 



dorapidities) 90°(?7 = 0), 40°(r/ = 1), and 4° (77 = 3.4). 
The analysis of the elliptic flow behavior together 
with the particle spectra makes it possible to deter- 
mine the extent to which the fall off of the integral v 2 
values can be attributed to a softening of the particle 
spectra going to forward rapidity. 

2 Experimental Details 

Elliptic flow measurements characterize the az- 
imuthal dependence of the particle production with 
respect to the reaction plane. For any given colli- 
sion, an empirical reaction plane can be established 
by measuring the particle yield observed in detec- 
tors mounted symmetrically about the beam axis. 
BRAHMS used several different detector systems for 
this purpose: Two rings of modestly segmented Si 
strip detectors (7 segments per strip) were configured 
with 42 active segments surrounding the beam axis in 
each ring. Another ring using the same type of Si strip 
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detector was configured with six segments around the 
axis, and seven segments along the axis. A scintillator 
tile ring with six panels was also arranged around the 
beam axis. Depending on the location of the beam 
vertex, for collisions within ±20cm of the nominal 
vertex location the Si and tile detectors covered pseu- 
dorapidity values from -3 to +1.5. At a somewhat 
larger rapidity (as —3), Cherenkov radiators mounted 
to phototubes in one of the experiment's Beam-Beam 
counter arrays^ were also used for the reaction-plane 



determination. The general procedure for determin- 
ing the reaction plane and the corresponding reaction- 
plane correction factor is described by Poskanzer and 

[81 

Voloshin J . The reaction plane corresponding to the 
second moment of the angular distribution is found 
with 

1 y- Wj sin (20,) 
2 2^fw i coa(2tf> i y { ' 

where the sum is over all detector elements with ge- 
ometric weighting factors Wi and azimuthal angle 4>i ■ 
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Figure 1: v 2 (pt,?y) for pions (both charges), kaons (both charges), and protons (filled symbols). Theoretical 
values based on the hydrodynamic calculations of Hirano et al. [9] arc shown by the open circles, squares, and 
triangles. The open crosses correspond to r/m 4 calculations using the AMPT model [10]. The smaller (blue in 
color) points in the pion panels are from a mixed-event analysis, as discussed in the text. 



The angular correlation with respect to the reac- 
tion plane of particles detected in the BRAHMS mid- 
rapidity and forward spectrometers determine the ob- 
served v 2 values, with 

< s Hcos(2[0-vI/ 2 ])). (2) 

Here the angular bracket denotes an average over all 
events of a given class, such as all pions within a given 
range of transverse momenta. Since the BRAHMS 
spectrometers are small acceptance devices, all par- 
ticles detected in the mid-rapidity spectrometer have 

0°, while those detected by the forward spectrom- 
eter have r>« 180°. 

The true v 2 value, which is based on the actual 
reaction plane, is found from the u| value using the 



reaction-plane correction factor R, with 

v 2 = v° bs /R. (3) 

The reaction-plane correction for any given ring (or 
ring combination) is determined using two additional 
rings that are far enough apart so as to avoid au- 
tocorrelations. Considering rings A, B, and C, the 
reaction-plane correction factor for ring A is found 
as: 

_ / (cq S (2[^-^b ])) {cos{2 [^-M,c])) 
V (cos(2[*f-* 2 c ])) ■ {V 

The reaction-plane correction factor was typically in 
the range 0.18 < R < 0.25. 

Figure [T] shows the resulting values (closed sym- 
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bols) of v 2 (pt) for pions, kaons, and protons at the in- 
dicated pseudorapidities, selecting events in the 10% 
- 50% centrality class. The behavior for the three dif- 
ferent particle species shows very little change with 
pseudorapidity. The smaller (blue in color) filled cir- 
cles shown in the pion panels are the results of a 
mixed event analysis, where the reaction-plane an- 
gle used for the analysis of a given identified pion is 
taken as that obtained from the previous event where 
a particle of any type is detected in the same spec- 
trometer. Because of the limited angular acceptance 
of the spectrometers, the deduced v 2 (p t ) values from 
the mixed-event analysis reflect the integral v 2 be- 
havior. The scatter of the points with p t mirrors the 
statistical uncertainty in the measurements. The fall 
off of the average behavior, as indicated by the dashed 
lines, is consistent with the known fall-off of the inte- 
gral v 2 values going to forward rapidity. 

We have also measured the p t -dependent parti- 
cle spectra for the channels where the differential el- 



liptic flow has been determined. The average value 
of transverse momentum for pions in the 10% - 50% 
centrality class at 90° is (475 ±60) MeV/c, falling to 
(380±45) MeV/c at 4°. These values are obtained by 
fitting a power-law dependence to the respective par- 
ticle spectra. The softening of the particle spectrum 
going to more forward rapidities can have a signif- 
icant effect on the integral v 2 values, as illustrated 
in Figure [D The insert shows an assumed form of 
v 2 (p t ) that is taken as being constant as a function 
of pseudorapidity. Folding this behavior with the ex- 
perimentally observed normalized particle spectra for 
pions at 90° and 4° results in the solid and dashed 
lines of the main figure panel, respectively. The softer 
spectrum at the more forward pseudorapidity leads to 
the weighted v 2 distribution peaking at a lower mean 
p t value, resulting in a 22% smaller integral v 2 value 
than a mid-rapidity. Two-thirds of the observed in- 
tegral v 2 change is then attributable to the softening 
of the particle spectrum. 
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Figure 2: Illustration of how the softening of the particle spectra going to forward angles affects the integral 
v 2 values. The insert shows the general behavior of the differential v 2 (p t ) values for pions, which for this 
illustration is assumed to be independent of pseudorapidity. Folding this distribution with the corresponding 
experimentally observed pion spectra at 4° and 90° results in the dashed and solid curves, with integral v 2 
values of 0.046 and 0.036, respectively. 



The small change in the differential v 2 signal going 
from mid- to forward rapidity suggests a longitudi- 
nally extended region for the medium produced in the 
collision. Rapidity dependent changes in the radial 
flow or other rescattering behavior of the hadronic 
stage might then account for most of the fall off of 
the integral v 2 signal going to forward rapidities. 

Hydrodynamic calculations are able to reproduce 
the observed behavior. This is seen by the open cir- 



cles in Figure [T] which show the results of the hy- 
brid hydrodynamic calculations of Hirano et al. (°\ 
and private communication) that include dissipative 
effects of the late hadronic expansion stage. Good 
agreement is found with the experimental values. 
The open crosses in Figure [T] show the results of the 
AMPT model for pseudorapidity rj — A. At this angle 
the string melting mechanism that allows for a good 
reproduction of mid-rapidity results has been turned 
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Figure 3: a) v 2 {pt) vs. p t for pions and protons at 77=0, 1, and 3. b) Constituent quark scaled values, as 
discussed in the text. The lines indicates the band containing most of the previous experimental results near 
mid-rapidity. 



off (^, and private communication). Although the 
comparison is subject to poor matching of pseudo- 
rapidity, there is a suggestion that the longitudinal 
extent of the produced medium may be greater than 
that assumed in current AMPT calculations. 

The mid-rapidity elliptic flow behavior for many 
different systems have been found to follow similar 
constituent scaling behavior when the v 2 value per 
valence quark, scaled by the inverse of the Glauber- 
model eccentricity, is plotted as a function of the 
mean transverse kinetic energy per valence quark ^— ' . 
ft Figure [3^ we show the v 2 vs. p t behavior of pi- 
ons and protons for the three measured angles on the 
same plot. In Figure [3b these same data are now 
shown in terms of their constituent-quark scaling be- 
havior. The two solid curves mark the approximate 
limits established by the previous measurements. The 
current data show a very similar trend, although with 
a tendency to fall somewhat below the previous re- 



sults. This may reflect a systematic error in the mea- 
surement of the reaction-plane correction factor or in 
the assumed Glauber eccentricity which is calculated 
for the trigger- weighted center of the centrality range. 
An error in either of these quantities will tend to have 
a similar scaling effect on all of our results. 

3 Conclusions 

BRAHMS has measured v 2 (p t ,ri) for pions, kaons, 
and protons at ^/s NN = 200 GcV for the Au+Au re- 
action. The results indicate a longitudinally extended 
region is produced where the eccentricity of the cre- 
ated medium and the corresponding pressure gradi- 
ents remain remarkably constant. Hydrodynamic cal- 
culations with final stage dissipation are found to be 
in excellent agreement with the measured differential 
v 2 values. The data are consistent with constituent 
quark scaling over the covered pseudorapidity range. 
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